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Abstract The stochastic partial differential equation proposed nearly three decades ago by Kardar, 
Parisi and Zhang (KPZ) continues to inspire, intrigue and confound its many admirers. Here, we i) pay 
debts to heroic predecessors, ii) highlight additional, experimentally relevant aspects of the recently 
solved 1+1 KPZ problem, iii) use an expanding substrates formalism to gain access to the 3d radial 
KPZ equation and, lastly, iv) examining extremal paths on disordered hierarchical lattices, set our 
gaze upon the fate of d=oo KPZ. Clearly, there remains ample unexplored territory within the realm 
of KPZ and, for the hearty, much work to be done, especially in higher dimensions, where numerical 
and renormalization group methods are providing a deeper understanding of this iconic equation. 
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1 In a Nutshell 

The history of physics has been punctuated at seminal moments by the appearance of certain fun¬ 
damental equations (and associated models) which have vigorously propelled the enterprise forward, 
serving as an explosively rich departure point, generating a myriad of alternative perspectives, creative 
insights, surprising connections and, given its sustained impregnability, often remain for many years a 
sacred object of fascination and obsession to its dedicated disciples. Recent, obvious suspects in this 
regard include the quantum mechanical Schrodinger equation (equally well, its flipside- Feynman’s 
path integral formulation), or the wonderfully elusive Navier-Stokes equation governing fluid mechan¬ 
ics, by which may be gleaned the scaling secrets of turbulent flow, dynamically encoded in the whirling 
eddies drawn by da Vinci centuries ago. Within the domain of equilibrium statistical physics, the 2d 
Ising Model, with its tour-de-force algebraic solution by Onsager, followed by combinatoric, graphical, 
Grassmannian, Monte Carlo, as well as full-blown field-theoretic, scaling, and renormalization group 
treatments, represents an extraordinary legacy that continues unabated to this very day. Arguably, a 
non-equilibrium statistical mechanical analogue to Ising/Onsager is the iconic equation Jl] proposed 
a generation ago by Kardar, Parisi, and Zhang (KPZ); it captures the statistical fluctuations of a 
kinetically-roughened scalar height field h(x,t ): 

d t h = vS7 2 h + ^A(V/r) 2 + ^/Dg, 
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subject to a i) simple laplacian smoothening mechanism, ii) tilt-dependent local growth velocity, and, 
finally, iii) spatiotemporally uncorrelated noise ry(x,t), which works in competition with the preceding 
deterministic pieces of this highly generic, but stubbornly resistant stochastic PDE. Here, v, A and 
D are coarse-grained, macroscopic phenomenological parameters, calculable for any given microscopic 
model via the technical machinery of KPZ scaling theory, elucidated by Krug, Mcakin and Halpin- 
Healy (KMHH) [2JS], built upon earlier finite-size scaling ideas [4], but predicated ultimately on the 
central role of the key scaling parameter 9 = H 1 /“A, where A = D/2v, and a divergent, parallel 
correlation length = (\/ r A\t) 1 / z . The 1 st KPZ Revolution, 1985-1999, was characterized by exact [5] 
ii 7 s 9 \wnm- renormalization group min, mode-coupling m, and numerical investigations 
whose primary focus was pinning down the dimension-dependent saturation- 
width and transient-regime scaling indices: a and /?=“. Galilean invariance, a key underlying symmetry 
of the KPZ equation was guarantor of the sacred, dimension -independent index identity a + z = 2. 
There was, even at this early time, a genuine appreciation of amplitude ratios Eld, as well as 
indications of underlying, universal KPZ PDFs mm- Complementary experimental work was scant, 
but provocative- see early, comprehensive reviews [28U291I30] . This epoch closed with the beginnings of 
Finnish investigations m of kinetically-roughened KPZ firelines, key mathematical papers HT?l.nTTTI . as 
well as refreshing nonperturbative [34] and conformally invariant [351 perspectives, the former inspiring 
numerical rebuttals [551I3T1I551I55] of a stubborn, battered suggestion [3D], revived nevertheless shortly 
thereafter [41], that 4+1 might be the upper critical dimension (UCD) of the KPZ problem. Recent 
numerics |42ll43ll44ll45] appear to have buried this idea, but from the analytical side [461147] . there’s 
lingering suspicion that something nontrivial is afoot near this particular dimension. 

The 2 nd KPZ Revolution, commenced in 2000 with the spectacular physical insights of Prahofer 
and Spohn [45] on polynuclear growth (PNG), and the creative mathematical efforts of Johansson [45] 
on the single-step (SS) model, both making explicit connections to the related problem of directed 
polymers in random media (DPRM). These works established that earlier numerically observed 1+1 
KPZ height PDFs were simply zero-mean, unit-variance versions of the universal Tracy-Widom (TW) 
limit distributions d, well-known from random matrix theory- RMT. Here, flat (radial) KPZ growth 
scenarios corresponding to Gaussian Orthogonal (Unitary) Ensembles; i.e., GOE (GUE), respectively. 
Indeed, TW-GUE was the very same PDF at the heart of the Ulam problem [5X1152] . a purely math¬ 
ematical matter focussed on the fluctuation statistics of the longest increasing subsequence (LIS) of 
a random permutation of the first n integers. Indeed, at this time, Okounkov (53] had elegantly es¬ 
tablished the LIS-RMT equivalence, a connection suggested empirically from the numerical work of 
Odlyzko and Rains [54] , building upon nearly prehistoric efforts of Baer and Brock [55] , which indicated 
the inevitable centrality of TW-GUE limit distribution in both instances. Prahofer-Spohn emphasized, 
as well, the importance of the Baik-Rains (BR) F 0 limit distribution [56], with no RMT manifestation, 
dictating temporal correlations of the KPZ stationary-state. A separate, critical paper by Prahofer 
and Spohn [52 introduced the notion of the Airy process, bringing spatial correlations into the RMT 
context, and tightening ties between KPZ physics and TW mathematical communities. During this pe¬ 
riod, flameless Finnish firefronts continued to burn [551I5M5U] . mode-coupling [5T] and nonperturbative 
mm approaches developed further, and a suggestive Dutch chemical vapor deposition experiment 
appeared [64], providing evidence of an asymmetric 2+1 KPZ height fluctuation PDF. For a recent 
example of the nice interplay of KPZ/DPRM statistical physics and TW/BR mathematics, as well as 
a brief account of the foundational Ulam-LIS numerics of Baer and Brock, see [55]; note, too- [551157] . 

A 3 rd KPZ Revolution- in fact, a veritable KPZ Renaissance EEMI- was inaugurated in 2010 
with the elegant, exact, and nearly complete solution of the 1+1 KPZ/DPRM problem. This included 
the full time-evolution of the universal PDFs to their asymptotic TW forms, managed by independent 
researchers using complementary driven lattice-gas mm and replica-theoretic DPRM approaches m 
E3j, in the former instance relying heavily upon recently gained insights into the weakly asymmetric 
simple exclusion process m- The initial advance was for the KPZ wedge geometry with its TW- 
GUE connection. Shortly thereafter, the transient dynamics of the 1+1 flat KPZ/TW-GOE class was 
finessed by Calabrese and Le Doussal m using Ising-style Pfafhan tricks within the DPRM context; 
for DPRM-hard wall interaction, relevant to TW-GSE, see m- Lastly, the stationary-state statistics 
of the 1+1 KPZ equation and its relation to the BR-F 0 limit distribution was elucidated via replicas 
by Imamura and Sasamoto HZ]; later, more rigor m- These extraordinary analytical works were 
matched, simultaneously, by the experimental efforts of Takeuchi and Sano on turbulent liquid crystal 
kinetic roughening phenomena. Reliant upon the machinery of KPZ scaling theory, they first observed 
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TW-GUE statistics [71], then TW-GOE [501151] and, finally, identified a key experimental signature 
of the KPZ stationary-state [821165] . Experimental manifestation of KPZ BR-Fo remains, however, a 
pressing matter. Importantly, the present KPZ Renaissance has inspired a renewed interest in the 2+1 
KPZ problem, as well as other dimension-dependent issues, such as the existence of a finite UCD, 
and using Wilsonian renormalization group methods, hopes of elucidating the full KPZ phase diagram 
[83] . Universality of the 2+1 KPZ Class and its associated limit distributions (i.e., higher-dimensional 
analogs of TW-GUE, TW-GOE and BR-F 0 ) have been well-characterized [5¥1I5^1I5H1[57| . along with 
universal spatial and temporal correlators [55115^] . supplemented by secondary KPZ “patch” PDFs j5SJ 
[90] , On the mathematical side, the transmission has shifted into hyperdrive; a rapidly accelerating 
moveable feast- e.g., 


2 An Homage to PS 

In Figure 1, we tip our hat to Prahofer and Spohn [481 . plotting up the flat, stationary, and curved 
height fluctuation statistics for the KPZ equation itself , a counterpoint to the PNG model of the 
original opus, stacked up against the TW-GOE, BR-F 0 , and TW-GUE limit distributions. To be 
clear, in the first instance, our 1+1 KPZ Euler integration, with A = 20, uniform stochastic noise, 
and time step St = 0.01, was obtained via a system size L=250000, with averaging done over 4000 
realizations (thus yielding 10 9 data in our statistical sample...), and integrated through 2000 time- 
steps. To make comparison to TW-GOE, these 1+1 flat KPZ Class height fluctuations are shifted, 
rescaled and swapped for the order-1 statistical Tracy-Widom variable Xi) according to Xi = (h — 
Voo t)/{rt)P, the parameter r = \AP A, trivially related to KMHH 6 above, and /3=| the exact KPZ 
exponent in this dimension. With longer runs, suitably averaged, we have extracted the asymptotic 
growth velocity n oo =0.195935 for the KPZ equation, as simulated. Note that, in contrast to PNG, 
where the phenomenological parameters A, A and are known exactly, pinning down these quantities 
represents the primary challenge, beyond simply aggregating the height fluctuation statistics. In fact, 
as done in experiment mm and occasionally in numerics, we have matched the variance of the 
transient regime statistics (i.e., our 1+1 KPZ Euler pdf) to TW-GOE, fitting the second moment 
to (xi)=0.63805, which fixes U^=0.23111. The resulting numerical data set, when viewed against 
the known TW-GOE limit distribution, focusses attention upon the skewness and kurtosis, as well 
as the finite-time correction to the universal mean. In particular, our extracted value for the 1+1 
KPZ skewness, s=(5h 3 ) c /(Sh 2 )^ 2 = 0.279, as compared with accepted result si=0.2935. As emphasized 
recently [65] , r thus determined, can then be employed in the service of a genuine comparison to both 
the radial TW-GUE and stationary-state BR-Fo limit distributions. In other words, the KPZ scaling 
parameters, while model-dependent, are not subclass dependent- in a given dimension, the same r and 
Uoo apply to flat, curved and stationary cases, as noted already in experiment [51]. Thus, here, we take 
the numerical short-cut, using TW-GOE to determine , privileging the greater challenges posed by 
BR-F 0 & TW-GUE. 

Hence, we have also included in Figure 1 the centered, rescaled fluctuations of the height increment, 
Ah = h(t 0 +At) — h(t 0 ), casting the associated stationary-state KPZ statistics in terms of the 0(1) Baik- 
Rains variable: xo = (Ah — VooAty/^At) 1 / 3 , integrating the KPZ equation to a late time t 0 = 5000, 
then investigating the temporal correlations at a later time t a + At, with 1 < At < 25 -C t 0 , at a given 
site, sampling over the entire system of size L= 10 4 , then averaging over 10 5 realizations, again yielding 
a statistical data set with 10 9 points. Here, for the skewness and kurtosis, we find 0.348 and 0.261, 
compared to Prahofer-Spohn’s oft-quoted BR-F 0 values: So=0-35941 and fco=0.28916, respectively, for 
these quantities. Additionally, our 1+1 KPZ Euler estimate here for the variance, 1.1527, compares 
well to the known Baik-Rains constant (Xo)=l-15039 and, as is evident from Figure 1, the 1+1 KPZ 
Euler stationary-state statistics are in fine agreement with the BR-F 0 limit distribution, proper. Note, 
by construction, this PDF has zero mean; i.e., vanishing first moment. As discovered by Takeuchi 
[52] , both experimentally for the 1+1 flat KPZ class turbulent LiqXtal system as well as the PNG 
model numerically, there is a dip in the distribution skewness as a function of the dimensionless 
parameter r = At/t a , as it interpolates between its TW-GOE (r 1) and BR-F 0 (r <C 1) values, 
Si=0.2935 and so=0-35941, respectively. The universality of the KPZ Baik-Rains limit distribution, 
Takeuchi’s skewness minimum, and other aspects of the transient to stationary-state crossover for the 
1+1 flat KPZ class, were addressed in a detailed, careful study [65j examining many diverse models. 
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Fig. 1 1+1 KPZ Euler: Numerical integrations, for stationary-state (xo), transient flat (xi), and radial (\ 2 ) 
geometries, stacked up against relevant universal limit distributions Baik-Rains, Tracy-Widom GOE, and GUE, 
respectively. Inset: Skewness minimum, stm~ 0.22, exhibited by KPZ equation as experimental si gna ture of 
stationary-state statistics, recently observed in KPZ turbulent liquid-crystal kinetic roughening work |B2]. Note, 
in particular, the asymptotic approach from below to TW-GOE skewness s(r 3>l)=si=0.2935. 


In Figure 1 insert, we strike a final note capturing, in a single long run, the skewness minimum in 
the KPZ equation itself. Here, we have L=10 4 , as before, but now t a =5 and 0.02< At <1000, and 
find s m i n = 0.218(4) for «5.0, near prior findings (65], which indicated s m i n =0.225±0.005. We 

mention, in particular, that this quantity, being the ratio of moments, requires no knowledge at all of 
the KPZ scaling parameter F, and emphasize its great value as an experimental signature of the KPZ 
stationary-state. 

Lastly, we include within Figure 1, a KPZ Euler integration relevant to the radial fluctuations in the 
TW-GUE class. We mention that, aside from some early works [1011110211103] . and more recently [104j . 
there has been little direct effort on numerical integration of 2d KPZ equation in polar coordinates; 
see, too [10511106(1107] . Indeed, all work on this subclass, aside from radial Eden model simulations 
[10811109] . have resorted to difficult, somewhat frustrating pt-pt simulations of various KPZ/DPRM 
models [1101165] in what is, effectively, constrained wedge geometries. The frustration arises because, in 
contrast to simulations for the flat KPZ subclass where all substrate points contribute to the ensemble 
average, the pt-pt Monte Carlo yields a few datum only per realization. On a practical level, this 
results in thousand-fold less payoff, given present-day CPU capabilities and equivalent run times; that 
is, probablilities down to 10” 6 vs 10 —9 , for TW-GUE vs. GOE, respectively. In any case, we report 
here results for KPZ Euler that rely neither on polar coords, nor a constrained pt-pt Monte Carlo, but 
rather upon an interesting numerical trick, built upon expanding substrates [59]. We refer the reader 
there for technical details, but the basic idea is to perform the simulation in the flat geometry, starting 
with a tiny substrate of system size L = Fo, but increasing the size of the substrate by 17 sites per unit 
time, stochastically duplicating column heights of randomly chosen points. This was, effectively, what 
happened in the polar coordinate KPZ integrations, in any case, since the kinetically-roughened frontier 
is expanding radially and its discrete mesh must be enhanced; i.e., its population of points grown in 
time. In practice, one choses Lq = 17 so, after an elapsed time t , this expanded substrate has grown to a 
size L = fit. In the original work, this method was used to check TW-GUE predictions for three specific 
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stochastic growth models: RSOS, SS and Etching, with typical values L 0 = Q = 20 and t = 500. A key 
ingredient involves the addition of a logarithmic time-dependence to the standard KPZ scaling ansatz: 
h = Voct + (rtfx + with C an independent stochastic quantity. As a warm-up, we have checked 
this methodology in the directed polymer context, using our well-characterized Gaussian g 5i DPRM 
from previous work |65| , where we performed a traditional, but quite labor-intensive zero-temperature 
pt-pt Monte Carlo simulation, with extremal paths of length t = 300, averaging over 10 s realizations 
of the random energy landscape, obtaining ((X2), (xi)> s 2, fo) = (—1-77097,0.8167,0.2304,0.0924), 
quite close to Bornemann’s state-of-the-art numbers mu. Here, we report an expanding substrates 
implementation of this same gaussian polymer model, using L o =20, t=2000, and averaging over 25000 
realizations, which generates 10 9 data at the final time slice, resulted in quite similar estimates (- 
1.77135,0.8311,0.2239,0.09373). In fact, these new g5\ DPRM results seem superior to KPZ growth 
model implementations on expanding substrates; specifically, in the admittedly difficult task of pinning 
down the kurtosis, which appears systematically underestimated, esp., the Etching model- see Ref|89|- 
Fig 2b, which indicates k? ~ 0.07-0.08. 

Performing the same trick for 1+1 flat KPZ Euler, however, presents its own difficulties since long 
simulation times are not easily accessible. Nevertheless, with Lo = 20 and 10 4 time-steps, which implies 
fully expanded substrates of size L=2xl0 5 , and averaging over 5000 realizations, we obtain here the 
following 1+1 KPZ Euler estimates for TW-GUE moments: (-1.77682,0.8108,0.2121,0.0826), and show, 
in Figure 1, the associated probability density in X 2 , measured against TW-GUE. We have taken the 
height fluctuation PDF at the final time slice (f=100), rescaled the width of the distribution with 
the model-dependent parameter 1^=0.23111, determined precisely by the earlier fit to TW-GOE, and 
then shifted the curve horizontally to the asymptotic, but model-specific, mean value (X2)=-K77682 
dictated by a full, finite-time scaling analysis. In any case, the agreement between our new 1+1 KPZ 
Euler integration and the known TW-GUE trace is reassuring. That said, we mention a slight residual 
underestimate of our 2d radial KPZ kurtosis /C 2 , seen previously for the deposition growth models. 


3 d=l+l KPZ Class Experiments: Stationary-State Metrics via WBC 

3.1 KPZ “Patch”-PDF I: Roughness Distributions for 1 /f 2 Noise &: Wiener Paths 

Beyond classic Tracy-Widom and Baik-Rains limit distributions capturing universal properties of ki- 
netically roughened edges in the 1+1 KPZ Class, there exist additional PDFs of relevance which 
characterize various intriguing aspects of the KPZ stationary-state statistics. We discuss in the follow¬ 
ing two such distributions which, experimentally and numerically, can be accessed via locally-averaged, 
so-called “patch” PDFs tied to sampling at small scales. These PDFs can play an important supporting 
role in the analysis of stochastic growth phenomena and have been conceived with the exigencies of 
the KPZ experimentalist well in mind, privileging window boundary conditions (WBC) over the theo¬ 
rist’s favorite- PBC. For an early pre-TW/BR indication of these practicalities, dating to the 1 st KPZ 
Revolution, we mention the efforts of Racz and Plischke mu. a longstanding statistical quantity of 
interest has been the squared-width: w 2 =(h 2 )-(h} 2 ; here, angular brackets signifying a spatial average 
over a tiny microscopic patch of linear dimension i, smaller than the finite dynamical correlation length 
£||. With an ensemble average over snapshots in the experimental system, the task is to aggregate the 
statistics of w = ( w 2 — (w 2 ))/cr^ 2 , which is the centered, rescaled, dimensionless order-one fluctuating 
statistical variable at the core of the universal roughness distribution- Pre,(lo). Recently, this metric 
was applied in 2+1 KPZ Class vapor deposition experiments involving organic [88j and semiconductor 
m thin films. In the former instance, Halpin-Healy and Palasantzas compared experimental data di¬ 
rectly to their numerical 2+1 KPZ Euler RD, using WBC “box”-sizes i -€£ ppz , for which the statistics 
were stationary, yielding nearly constant quantities. The 2+1 KPZ Euler RD was characterized by a 
skewness s=2.03 and kurtosis fc=7.11, while the organic thin-film experiment yielded measured val¬ 
ues: (s, fc)=(2.12,7.77). Since the roughness distribution, by construction, is a zero-mean, unit-variance 
PDF, s & k become the two essential quantities to pin down. Thus, the agreement is quite suggestive 
and given the great demands of thin-film experiment, as well as the relative paucity of temporal data 
there, reveals Prd(oj) to be a welcome, supplementary measure jll2l of 2+1 KPZ kinetic roughening. 

Motivated by these 2+1 KPZ vapor deposition experiments, we have revisited the 1+1 KPZ Class 
turbulent liquid-crystal work of Takeuchi and Sano MU, examining the data from the perspective 
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Fig. 2 Experiment, Numerics & Theory: 1+1 KPZ Class RD, semilog plot. Comparison- 1+1 KPZ Euler 
numerical integration, Gaussian DPRM, turbulent liquid-crystal data \HMMH and the roughness distribution 
of Antal et al. for 1 /f 2 noise [115] . Insets: a) Linear plot. Near the peak, thin KPZ Euler traces move upwards 
for decreasing box-size, b) RD skewness and kurtosis, as a function of patch-size i. Dashed horizontal lines 
indicate associated AntaEw values. In experimental system, £ is measured in pixels, corresponding to 4.74 
microns. 


of patch statistics and the local roughness distributions. Conveniently, our efforts here have been 
preceded by two extremely helpful theoretical papers. The first, by Foltin et al. ms, calculates the 
exact roughness distribution of lcl random-walk Edwards-Wilkinson (EW) interfaces 11141 . assuming 
PBC. We remind the reader that d=l+l is a very special case and, as regards the steady-state statistics, 
the KPZ nonlinearity is irrelevant in this dimension, with both EW and nonlinear KPZ proper sharing 
the same stationary-state measure. While this statement is typically voiced with the thermodynamic 
limit in mind, it is also true for a finite system with PBC; see Fokker-Planck discussion, next section. 
In fact, these authors numerically verify the 1+1 EW/KPZ steady-state equivalence in this PBC 
context, extracting the roughness distribution of catalytic reaction model, honorary member of the 
1+1 KPZ Class; see m -Fig 4 there, which reveals very fine KPZ agreement (in fact, better than 
EW!) with their exact RD for periodic random walks. That’s great, but we consider patch-PDFs using 
experimentally-dictated WBC. Serendipitously, later work by Antal and coworkers m addresses 
precisely this case; i.e., finite Brownian paths with WBC. Truth be told, their paper considers the 
much more general issue of RDs of periodic and nonperiodic 1 //“ noise. Interestingly, their results 
establish that, generically, the respective RD PDFs are distinct, PBC vs. WBC- see [115j -Fig 6, which 
makes this quite clear. Here, a =2 corresponds to 1+1 random-walk EW interfaces, and they obtain the 
exact roughness distribution for WBC; see Ref. [115j -sect. IV.A-B & Figs. 4 & 5, which illustrates the 
relevant Brownian RDs. In Figure 2 above, we compare this Antal 2 W RD with Takeuchi-Sano liquid- 
crystal data sets, as well as our newly-obtained 1+1 KPZ Euler WBC roughness distribution- here, 
A=20, with system size L=10 5 , and averaging done over 1000 runs. The main plot is semi-logarithmic, 
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fully revealing the long tail, with probabilities down to 10~ 5 ; inset left-Fig. 2a, illustrates the straight- 
up RD with linear scales. Additionally, statistics for a Gaussian polymer model, g5\ DPRM discussed 
previously, with normally distributed site energies has been included to indicate universality of our 
1+1 KPZ/DPRM results. Finally, Fig. 2b records the ^-dependence of the 1+1 KPZ RD skewness and 
kurtosis which, with diminishing box-size em, is observed to approach values characteristic of the 
Antal 2 m RD, calculated by us to be (s, fc)=(2.55,10.27), both s and k distinctly larger than 2+1 KPZ 
RD counterparts [5Hj: (2.03,7.11). On the experimental side, the LiqXtal RD skewness rides slightly 
above 1+1 KPZ, but the kurtosis suffers a bit more from the noisy statistics at the end of that long 
right-tail, which we extend nearly 12 standard deviations. We mention, as an aside, digging that far 
into the right tail of AntaU™ RD trace demanded five terms in the Laplace transform series, rather 
than the three of the original paper m ■ In this regard, we see pure exponential behavior for WBC, 
as is already known for periodic boundary conditions, though with a smaller prefactor than the n 2 /6 
characteristic of Antal 2p , appropriate to PBC [113 ]. All this, of course, contrasts with the 2+1 KPZ 
Euler roughness distribution [88] which, for WBC, appears to possess a stretched exponential tail with 
exponent 7/8. 


3.2 KPZ “Patch”-PDF II: Extremal Statistics & Majumdar-Comtet Distribution 

Our second patch-PDF focusses on the extreme-value (EV) statistics of the height fluctuations, provid¬ 
ing a kinetic roughening context in which to study extremal behaviors of correlated random variables. 
In this setting, a natural statistical variable is the maximum-relative-height (MRH), m = h m ax ~ h, 
measured with respect to the mean interface position in the patch. Extremal MRH fluctuations on the 
global scale, in large systems assuming PBC, was first examined numerically m for the Edwards- 
Wilkinson case, but then solved exactly in 1+1 dimensions in beautiful work by Majumdar and Comtet 
[11811119] , Somewhat surprisingly, they discovered that the Airy PDF, famously implicated in dictating 
areal statistics beneath Brownian bridges, was at work here as well, controlling MRH fluctuations of Id 
random-walk EW interfaces. As these authors emphasize, and check, the Airy distribution also governs 
MRH statistics for the 1+1 KPZ stationary-state [1181111911120) . This is not unexpected, since the KPZ 
nonlinearity generates boundary terms which vanish, yes, in the asymptotic limit, but also in finite 
systems with PBC; in other words, for PBC, the same stationary measure solves the Fokker-Planck 
equation for both 1+1 EW and KPZ, a well-known special feature unique to this dimensionality. Un¬ 
fortunately, as pointed out by Majumdar and Comtet, this bit of magic does not occur for the WBC 
relevant to us. Additionally, they calculate a quite distinct PDF relevant for 1+1 EW interfaces with 
“free” (i.e., Neumann) BC, which they christen “F-Airy”; see, Ref. [118] -Fig 1, for a comparison of 
these two universal PDFs. 

Given the missing magic, as well as the additional potential wrinkle of Neumann versus strict 
window boundary conditions, there is no reason necessarily to expect that this new F-Airy distribution 
will characterize the 1+1 KPZ Class WBC situation too- for key details, we refer to Ref. m -sect. 
5. Interestingly, despite these caveats, our analysis suggests that the Majumdar-Comtet PDF does, 
in fact, capture the stationary-state MRH fluctuations of 1+1 KPZ Class, thereby preserving the 
KPZ-EW steady-state equivalence and Id random-walk interpretation, WBC notwithstanding. We 
direct the reader’s attention to Figure 3 where, employing WBC, we examine MRH statistics in local 
patches of size £ <g£p pz showing turbulent LiqXtal data sets, and making explicit comparison to 1+1 
KPZ Euler integration and g 5i DPRM simulations analyzed from this perspective. Again, the main 
point is that by restricting ourselves to sufficiently small £ values, we are actually accessing the 1+1 
KPZ stationary-state. We include, too, the Majumdar-Comtet F-Airy distribution, proper. All this 
is presented in a zero-mean, unit-variance rendering to highlight the skewness and kurtosis of these 
PDFs. The main plot is semi-logarithmic, revealing the F-Airy tails in their full glory; for large, positive 
p = (m — (m))/cr m , the behavior is known [118] to be quadratic. That is, a gaussian tail, seen in the 
central inset, which uses linear scales. A second inset, upper right, reveals the ^-dependence of the 
skewness and kurtosis of our 1+1 KPZ Euler MRH patch-PDFs. Here, dashed lines correspond to the 
values (s, fc)=(l.ll,1.69) characteristic of the Majumdar-Comtet distribution, which we have extracted 
via its 10-term confluent liypergeometric representation [119] , Our thick KPZ Euler curve corresponds 
to £=56, where we have (s, A;)=(1.01,1.48). As the patch-size shrinks, however, greater stationarity is 
insured and agreement improves- see thin traces, corresponding to £=48,32,16; this is especially evident 





















in the center inset, near the peak of the F-Airy distribution. Hence, for £=32, our KPZ Euler MRH 
patch-PDF possesses (s, fc)=(1.05,1.57), which increases further to (1.07,1.62) for £=16, clearly heading 
towards calculated F-Airy values. We find these vanishing differences rather suggestive and note, in 
passing, that the Airy PDF appropriate to 1+1 KPZ/EW with periodic boundary conditions [1211)1221 
has (s,fc)=(0.701,0.560), certainly quite distinct from Majumdar-Comtet F-Airy, relevant for WBC, 
which we illustrate in Figure 3 and see in the LiqXtal experiment. 

Of course, the MRH distribution governing stationary, extremal height fluctuations of 2+1 KPZ 
surfaces presents a substantial theoretical challenge, as there are no ties at all to an EW counterpart 
11731 in this higher dimension, regardless of BCs. Nevertheless, there has been recent progress on 2+1 
KPZ Class extremal-height numerics, where explicit comparison has been made with MRH patch- 
PDFs obtained from thin-film stochastic growth experiments IHHUSDI- An interesting question there 
then concerns the higher-dimensional analog of the 1+1 KPZ Class gaussian F-Airy tail; at first 
glance, the 2+1 KPZ Euler MRHD [88] with (s, fc)=(0.884,1.20), though surely different, bears some 
familial resemblance to the Majumdar-Comtet distribution. Even so, it would be rather unexpected 
were the 2+1 KPZ MRH right-tail asymptotics also dictated by a precise gaussian fall-off. In any case, 
this would be quite a difficult matter to resolve via numerics. 


4 Universal Limit Distribution: 3d Radial KPZ Class 

In previous works [84U85] , attention was focussed on large-scale numerical simulations of multiple mod¬ 
els (RSOS, EDEN, driven-dimer, DPRM), as well as a direct Euler integration of 2+1 KPZ equation 
itself, establishing universality of the 2+1 flat KPZ Class, manifest in a limit distribution P^-p Z (£) 
possessing universal moments ((£i), (£i ),s,k)= (-0.85,0.235,0.424,0.346), with £ = (h — Voot)/(dt)K 
We draw particular attention to [55]-Section III, Figs. 6 &7, where a refined portrait of this character- 
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Fig. 4 Universal Limit Distribution: 3d radial KPZ Class, via expanding substrates. Insets reveal finite-time 
scaling analyses necessary to elicit asymptotic estimates for universal moments. Dashed horizontal lines indicate 
known values extracted from independent 3d pt-pt DPRM, SHE, and Eden simulations Eng. 


istic PDF, higher-dimensional analog of TW-GOE, is presented. See, too, (Mj-Table I & Fig 2 there, 
where an exhaustive analysis, employing KPZ scaling theory [2], built with the KM toolbox, permit¬ 
ted determination of all relevant model-dependent parameters. In particular, for the highly nonlinear, 
A = 20, 2+1 KPZ Euler integration, the asymptotic growth velocity = 0.17606, and key scaling 
parameter 9 = A 1 t a \ = 1.192xl0 -3 . Indeed, with knowledge of these parameters, it was also possible 
to isolate the limit distribution dictating the late-time 2+1 KPZ stationary-state statistics, analog of 
1+1 KPZ Class BR-Fo- This was done for 2+1 KPZ Euler, gDPRM, and RSOS stochastic growth 
models, demanding quite long integration times, but resulted in a convincing, composite portrait. The 
KPZ stationary-state is quite important as a subclass because it represents the natural contact point 
with field-theoretic, renormalization group and mode-coupling analyses; here, Halpin-Healy extracted 
the universal variance of the 2+1 KPZ stationary-state statistics, obtaining (£^) = 0.464, analog of 
the 1+1 Baik-Rains constant (Xo) = 1.15039. When properly translated [851 into the context of a 
recent field-theoretic, Wilsonian RG calculation of Kloss, Canet and Wschebor [55] , this becomes the 
universal renormalization group amplitude i? 2 +i = 0.944 ± 0.031, in fine agreement with the KCW 
value 0.940. Hence, we see that (^) = 0.464, where field-theoretic RG meets Monte-Carlo numerics 
has become, suddenly, the best measured quantity of the 2+1 KPZ problem. We remain hopeful that 
the field theorists will someday succeed in crafting their own renormalization-group portrait of the 
entire limit distribution, or perhaps the skewness therein, known to be s 0 = 0.244; see Ref [55] , Sect. 
IV and Fig. 9 in that paper. Interestingly, the 2+1 KPZ stationary-state has also permitted a solid 
multi-model (gDPRM,RSOS,KPZ Euler) estimate /32 +i= 0.241(1), in line with prior gold-medal studies 
VMFM at the 3-digit precision level. 

A more stubborn challenge involves unearthing the higher-dimensional KPZ analog of TW-GUE; 
i.e., the limit distribution characteristic of the 3d radial KPZ problem. The 2010 tour-de-force exact 
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solutions notwithstanding, solving the KPZ equation in 3d spherical coordinates, or the related 3d pt- 
pt stochastic heat equation (SHE) is a humbling prospect, indeed. Fortunately, recent numerics have 
helped chart the course, framing the essential features and providing a target for future analytical 
efforts. Halpin-Healy |Mj simulated 3d extremal paths connecting far corners of a cube filled with 
exponentially distributed site energies; the underlying pt-pt energy fluctuation limit distribution was 
measured to possess universal variance, skewness and kurtosis: 0.291, 0.331, and 0.212, respectively. 
Later authors [SB], studying three distinct 3d pt-pt KPZ growth models, among them single-step 
(SSC) and on-lattice version of Takeuchi Eden D, find values in the range S 2 « 0.32-0.34 and k 2 « 
0.20-0.22, providing additional evidence for these characteristic, universal quantities. Tucked away in 
his doctoral dissertation m- Table 7.2, Prahofer presciently records S2=0.323(5) and fc 2 = 0 . 21 ( 4 ) for 
his own early 3d pt-pt PNG simulation. Subsequent work iSS], see particularly sect. II: Table I and Fig. 
4 there, which involved direct numerical integrations of the constrained 3d pt-pt SHE, as well as three 
additional 3d pt-pt DPRM models, has established ((£ 2 ), (£|)> s 2 j k 2 )= (-2.28,0.325,0.328,0.214) in this 
dimension. Here, we perform a 2+1 KPZ Euler integration on expanding substrates with 12 = 10, 
complementing our earlier efforts on the 3d pt-pt SHE with multiplicative noise. Results are indicated 
here in Figure 4, with insets demonstrating that the variance, found to be 0.326, is dead-on, the 
skewness and kurtosis well within accepted values while subleading corrections render 

extraction of the universal mean a bit more challenging. Nevertheless, successive fit lines monotonically 
approach a well-defined boundary (i.e., envelope), with intercepts converging to (£ 2 ) ~-2.32, a quite 
decent value. 


5 Diamonds in the Rough: The Hierarchical <C>DPRM Revisited 

5.1 From Gauss to Gumbel back again 

The dimensional behavior of the Euclidean KPZ/DPRM problem can be found, in nearly its full glory, 
within a directed polymer implementation on b— branched, diamond hierarchical (•()) lattices; see Figure 
5, lower left, for geometric reminder. At the heart of the <£>DPRM one finds a simple recursion relation 
[12511126] . whose fixed point function, see Fig. 5 for examples, is the much-desired, sought-after limit 
distribution: 


pOO 

■ poo 

/ Pn+l = 

/ Pn* Pn 

k 

Uz J 


Here, P n represents the <)>DPRM bond-energy PDF at the n th generation, and the asterix indicates 
convolution. The pioneering effort in this regard, dating to the 1 st KPZ Revolution, is due to Derrida 
and Griffiths [125] , an elegant work which quickly inspired an analysis of its many-dimensional gener¬ 
alization [ 126] , addressing the question of a finite UCD in the KPZ context. Indeed, it was upon these 
lattices that universal DPRM PDFs were first examined- see, esp., Halpin-Healy jTD]-Fig. 8 there, which 
shows the asymptotic evolution of underlying distributions, as well as Derrida [127] - Fig. 3, providing 
a representative snapshot of a, now, very familiar (or, at least very familiar+oofcmg....) asymmetric 
curve. Following up on the <0>DPRM, Roux et al. [128] proposed an index relation r/ = 2r /, connecting 
exponents dictating the stretched exponential tails of the universal DPRM energy PDF. With Zhang’s 
saddle-point formula for the left-tail exponent, ry=l/(l-w) and the known 1+1 KPZ/DPRM result 
Pkpz = udprm = 5 in hand, this led to the conclusion, addressed numerically j5B] back in 1991, that 
77 = | for the long, TW Airy tail, but also predicted t)'=3 for the short, computationally intractable, 
steep opposing tail. These values, of course, were understood later by Prahofer and Spohn, as char¬ 
acteristic of the underlying universal Tracy-Widom limit distributions. Meanwhile, the <)DPRM lay 
dormant [indeed, for some time....] until Monthus and Garel (MG) returned to this rich, deceptively 
simple toy model. In an intriguing paper [129 ]. they i) announced a revised tail exponent conjecture: 

77 ' = d eS ri 

for the DPRM problem, ii) provided a saddle-point analysis tailored to and specifically valid for the 
<>DPRM, where d e g = 1 + log 2 & serves as the fractal dimension for 6 -branched hierarchical lattices, 
but also, iii) proposed a simple back-of-tlie-envelope argument extending the result to the standard 
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Fig. 5 Universal Limit Distributions: ^>DPRM Class- Lattice branching parameters 6=1, 2.82826, and oo, 
the first and last values yielding normal Gaussian and standard Gumbel PDFs, respectively. Inset, upper left: 
6-dependence of the KPZ scaling exponent lower left: Geometric gnomon, 6-branched hierarchical lattice. 


Euclidean DPRM, where d e g = d+ 1, literally. Lastly, the curious researcher might also compare MG- 
Fig 3a to |26j-Fig 9; the similarity is not incidental, but rather intentional and, indeed, quite germane 
to the discussion that follows. 

As for the Monthus-Garel conjecture, we direct the dedicated reader to Sections 2, 5.1, as well 
as Appendix A of this interesting work [129] . which contrasts their own rare events analysis to the 
earlier argument of Zhang [29] that had first fixed the TW Airy tail exponent ? 7 =l/(l-w). We also draw 
specific attention to MG’s numerical confirmation of Zhang’s formula for the hierarchical <£>DPRM; 
esp., MG-Table 1, where they explore numerous integer values of the branching parameter, running 
from 6=2, where w=0.299, so 77 = 0 .143, to 6=24, where w=0.123 yields 77 = 1 .14. Unfortunately, the more 
challenging task of actually testing the MG conjecture, itself, was sabotaged by the large corrections 
to scaling attendant to rj. Nevertheless, the MG data appear “compatible” with their suggested tail 
exponent identity. Finally, the authors present an impressive group portrait of the associated limit 
distributions for 6 in this range; see, MG-Fig. 3b. 

We have revisited the 6 -dependence of the <C>DPRM universal PDF and, show in Figure 5, our own 
portrait focussing on three seminal values of 6 . That is, the limit 6 —* 1, i.e., d= 0, the model degenerates 
to simple sums of uncorrelated random variables, the variance grows linearly, w=l/2, and the CLT 
demands symmetric Gaussian; in other words, r]=r] l =2, and s=k= 0. Conversely, as 6 —> 00 , where there 
are an infinite number of branches, one effectively selects the extremum of an asymptotically large 
set of iid’s. Indeed, in this limit, one retrieves the textbook, extremal statistics asymmetric Gumbel 
distribution [130] . which possesses a skewness SG= 12 -\/ 6 £( 3 )/ 7 r 3 , involving Apery’s constant, an excess 
kurtosis kc=^§-, a simple exponential for the long, left-tail, and a deep-diving, double exponential for 
the steep right-tail. Thus, the ^DPRM provides a convenient continuous model interpolating between 
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Fig. 6 Hierarchical (•) vs. Euclidean (•) DPR.M: sfc-plot. Power-law fit to the final dozen points of the 
<C>DPRM data set, indicated by the colored trace ( ), nearly pierces the bull’s-eye at the upper right corner 

of the plot frame, its expected destination, where reside the standard Gumbel values ( sg , /cg)=(^^C(3), ^). 
Inset- Upper left: b —dependence of the <£>DPRM skewness and kurtosis. They cross, for 6«2.82826, at the 
common value s=k~0. 49375, indicated by the dashed line. Lower right: Running exponent p e g in the <C>DPRM 
power-law fcs p -relation. 


normal Gaussian (e - ^ ) and standard Gumbel (e^ _e? ) limit distribution behaviors for b £ [l,oo), as 
we emphasize in the figure. The key point here is that, as b grows, the long left-tail swings up and 
•v>DPRM exponent 77 drops from 2 to 1, evolving from quadratic to linear, while the index rj starts 
at 2 but then diverges without bound, transforming the steep right-tail into a double exponential 
e -e in the infinite -6 limit. A final inset to Figure 5 reveals the 6 -dependence of the characteristic 
KPZ scaling exponent w( 6 ), which controls the growth, at the n^-generation of the lattice, of the 
ODPRM PDF width S n ~ 2 nul , as well as the vanishing deviation of the distribution mean 7 n , from 
its asymptotic value; the latter behavior not shown but, interestingly, fully consistent with Euclidean 
KPZ expectations [131 . 

Here, we concentrate specifically on the <£>DPRM skewness and kurtosis; see Figure 6 , where we 
cast s vs. k, leaving 6 momentarily behind the scenes, crafting a parametric plot. Presenting the data as 
such, stripped of its hierarchical lattice origins, will permit us salient comparison to Euclidean DPRM 
results later. Even so, within the figure inset, upper left, the explicit dependences, s( 6 ) and k(b) are 
indicated, where we have chosen a logarithmic scale for the abscissa, in accordance with expectations 
regarding d e s(b) for the <0>DPRM. Hence, this inset reveals the behaviors of s and k as a function of 
the branching parameter 6 ; i.e., effectively, dimensionality. In particular, while s and k both vanish in 
the Gaussian limit (6 = 1), the skewness rises more rapidly than the kurtosis initially, scaling roughly 
as the square root. By contrast, the kurtosis appears to have positive curvature at the inset origin, but 
quickly straightens out to something essentially linear, though with slightly negative curvature beyond 
the curious (inflection?) point where the kurtosis and skewness intersect, sharing the common value 
s*=k*~0 .49375 at 6«2.82826; recall, again, Fig 5. Of course, negative curvature in k is expected, given 
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the anticipated infinite 6 limit; i.e., standard Gumbel, with s< 3 =1.14 and /cg= 2.4. We emphasize that 
our measured s and k values are asymptotic, invoking the characteristic KPZ finite-time correction, 
which vanishes as t~ 2ul . Returning to Figure 6 proper, one senses the possibility of a power-law relation, 
k ~ s p , between the <0>DPRM skewness and kurtosis; perhaps, parabolic. This is investigated in a second 
inset, placed lower right, where we plot the running exponent p e s=dlnk/dins. Indeed, for increasing 
“dimensionality” (i.e., b f), as both s and k grow beyond the common intersection s*=k* « | , heading 
towards their respective standard Gumbel values, the index p k 2~. By contrast, in the opposite limit, 
as b heads towards unity, where the limit distribution becomes Gaussian with vanishing skewness and 
kurtosis, we see a different, quite distinct exponent p ss 1 + . That said, it is manifest already that at 
6=2, where s=0.368, the index p ss 1.79, well on its way towards quadratic behavior. 

Finally, Figure 6 proper also includes individual data points associated with our measured (s, k) 
values for the d+1 Euclidean KPZ/DPRM problem, in transverse dimensionalities d= 0-6. We note, in 
particular, that for unconstrained 3+1 DPR.M extremal trajectories, with Gaussian distributed random 
energies on the sites [85] , we find s=fc=0.51 _ ; see discussion below. Thus, somewhat serendipitously, 
it appears that 3+1 KPZ provides that instance in which the skewness and kurtosis slide past each 
other numerically; in fact, within the error bars, it is the same such value exhibited by the 6=2.9 
ODPRM. Of course, the coincidence of s and k hasn’t the slightest mathematical significance. It is, 
rather, just an idiosyncratic feature of the 3+1 flat KPZ Class; another being the existence of a finite- 
temperature phase transition separating entropic (T > T c ) and extremal ( T < T c ) DPRM wandering 
m- Nevertheless, that the coincident value is nearly the same, s*=k* ~ \. for both the Euclidean and 
hierarchical DPRM is, quite frankly, a little surprising. Additionally, it raises many natural questions 
regarding the fundamental characteristics shared (and not...) of these similar, yet topologically distinct, 
extremal path lattice models. 


5.2 The Many-Dimensional DPRM Fate of d—o o KPZ 

We comment further here regarding Figure 6 and associated results, as well as the higher-dimensional 
KPZ context within which they sit- 

• Prior studies |125lll26U127U128Ul29| of the <0>DPRM focused, initially at least, upon the exponent 
w and had confirmed (one might say, in part, lamented...) the slightly “low” value w=0.299 found 
for b=df s =2, given the exact result w=|, well-known for the 1+1 KPZ/DPRM problem. A glance 
at Figure 5 insert, however, reveals that there is, of course, a value of the branching parameter for 
which this famous exponent is actually retrieved. Indeed, for 6=1.69, we find w=0.3338 _ AND, quite 
surprisingly (this is the main point....), s=0.2947 and fc=0.1956; i.e., values nearly coincident with TW- 
GOE results, Si=0.2935 and /ci=0.1652, respectively. An additional pay-off comes when we examine 
the asymptotic scaling behavior of the first and second moments, extracting the <0>DPRM analog of 
the TW-GOE universal ratio (xi)/(Xi ) 1,/2 = 0.9515, which marks the mean-to-width characteristic of 
the underlying limit distribution. Here, we discover 0.86. In other words, the 6=1.69 </>DPRM PDF 
provides a reasonably good match to its 1+1 KPZ/TW-GOE cousin. Since u is dead-on, the long TW 
Airy tail (rj = §) is captured, along with the distribution’s skewness and nearly, but not quite, its 
kurtosis. Agreed, the steep tail exponent, rf = d e sr] ss 2.64, is a bit shy of the Tracy-Widom value 
3 but, we note, with some irony, very much in line with existing heroic |261I133| . though admittedly 
insuffucient, numerical estimates of this stubbornly difficult, intractable quantity. Hence, in a zero- 
mean, unit-variance rendering, the 6=1.69 <0>DPRM PDF would sit comfortably amidst standard 1+1 
KPZ class model simulation data (e.g., RSOS, PNG, etc.) unless one plunged, at substantial numerical 
cost, deep into the tails. Nevertheless, as we will discuss in detail momentarily, the clear numerical 
distinction between these k values (esp., the fact that /+=i.69 > fcrw— GOEj despite the coincidence of 
the skewness s...) will be of great consequence. 

• Recalling recent 2+1 DPRM results [84], which have isolated the higher-dimensional analog of 
KPZ/TW-GOE, measuring s=0.424 and A;=0.346 for the generic pt-plane case, and pinned down the 
key index W2+i=0.241 [20111241 via a multi-model study [SO] of 2+1 KPZ stationary-state statistics, we 
notice quite similar values for the 6=2.3 //DPRM, where the exponent, 0.274, might be a little high, 
but the skewness, 0.423, and kurtosis, 0.372, certainly close to the mark. Given our findings above for 
the <(>DPRM sk relation, one cannot resist looking at Euclidean KPZ from this same vantage point. 
Indeed, naively fitting known 1+1 TW-GOE and 2+1 KPZ/DPRM [SO] values for the skewness and 
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kurtosis to a simple power-law yields fc=1.94s 201 ; i.e., a quadratic dependence. We encourage the en¬ 
gaged reader to dust off their calculator, check these numbers, and then “compute” k for TW-GUE 
s=0.2241. This is at least a little surprising, no? Next, consider the 3d pt-line DPRM [SSJ-Appendix, 
for which (s, fc)=(0.408,0.322), as well as the earlier discussed 3d pt-pt DPRM/KPZ problem, where 
(s,fc)=(0.328,0.214). In these instances, as well, the ratio k/s 2 =2~\ in fact, eerily so. Data points as¬ 
sociated with these additional pt-pt and pt-line KPZ subclasses have thus been included in Figure 6 , 
filling out the curve a bit. 

•Curious about these matters, we have gone on to higher dimensional hypercubic DPRM, measuring 
for the 3+1 flat case: w=0.1868, s=0.508 and fc=0.509, with system sizes Z/=500 and «10 n data. 
For a distinct, but related polymer with random site energies drawn uniformly from the unit in¬ 
terval, so-called u5 DPRM [84], we find quite similar values s=&;=0.515. As disclosed earlier, our 
main numerical observation regarding the 3+1 DPRM/KPZ problem is that already in this dimen¬ 
sionality, the limit distribution kurtosis has risen, numerically, to the level of the skewness. More 
to the point, however, we now supplement that arithmetic fact with our sfc-relation “prediction” 
s*=k*= 1.94 " 1 « 0.51. For the 4+1 Gaussian DPRM, we find w=0.152+0.004, s=0.577 and k= 0.688, 
respectively, for L=125. Things become progressively more challenging in higher dimensions, but for 
the 5+1 DPRM: (s, fc)=(0.596,0.73), and 6+1: (s, fc)=(0.623,0.805). These findings, when cast in terms 
of an sk-plot, suggest that the (fDPRM is quite distinct from its Euclidean cousin. We have illustrated 
this explicitly for the 6=1.69 hierarchical lattice, whose limit distribution possesses w=|, shares the 
TW-GOE skewness, but has a kurtosis, k=0. 1956(1), definitively larger than the Tracy-Widom value 
fci=0.1652. In other words, for small skewness, the sfc-curve of the <}DPRM rides well above the Eu¬ 
clidean DPRM data; however, the respective sk- traces eventually cross near the transcendent value 
s*=k* ss Since it is known that asymptotically, as b —> oo, the ODPRM limit distribution converges 
to standard Gurnbel- i.e., the point ( SG,kc ), upper right comer Figure 6 , proper- this suggests that 
d= oo KPZ, in its DPRM manifestation, is NOT simply Gurnbel. The sfc-trace of the many-dimensional 
DPRM is, at the start, distinct, but then rises precipitously at a different rate. Indeed, this sugges¬ 
tion is consistent with independent work [134] investigating directed polymers on Cayley trees [135] , 
effectively an infinite dimensional implementation of the KPZ problem. 

• Parting gesture: On a rather different note- Folding together in one hand the KK conjecture m , 
w = l/(d + 2), for the d+1 DPRM/KPZ problem, Zhang’s formula for the TW Airy tail, r] = 1/(1 —w), 
along with the MG tail-exponent relation leads to the hopelessly seductive, simple result that rj = 
(d + 2)/(d+l) and 77/ = d + 2 for Euclidean lattices; in other words, (77,77/) = (|, 3) known rigorously 
for the 1+1 KPZ/DPRM TW limit distribution, where d= 1, while for the 2+1 DPRM: (77,?;/) = (|,4), 
and in 3+1: (77,77O3+1 = (§ , 5) , etc. This is almost certainly not true, but frankly, nearly irrefutable 
numerically [133] (even admitting the loss of KK....) and, we must confess, an irresistible, if only 
convenient, rule-of-thumb. 


6 The Future 

We close with a passing observation, then a rallying call. Firstly, for a strict quadratic dependence, 
whereby p = 2, the quantity k/s p —> (66 4 ) c (Sh 2 ) c /(Sh 3 ) 2 , essentially the Derrida-Appert ratio [136 . 
From this perspective, our work here reveals that the Derrida-Appert ratio is, despite quite a broad 
range of underlying s and k values, nearly constant for the entire Euclidean KPZ problem, across many 
dimensionalities, as well as geometric subclasses: pt-pt, pt-line, pt-plane, etc. Indeed, the variation 
of this quantity is quite modest, running from TW-GUE: 1.86 to 6+1 DPRM: 2.07. Clearly, p / 2 
exactly. Nevertheless, we wonder whether a tightly-knit family of limit distributions may someday be 
implicated in these numerical conspiracies, or otherwise explain the findings of Figure 6. Until then, we 
encourage mathematically-inclined colleagues to revisit the <)DPRM and related matters. Any rigorous 
statement (e.g., a bound) regarding the tail exponents of its 6—parametrized limit distributions would 
be a great boon, certainly. Of course, an exact, full solution for a single 6-value wouldn’t be turned 
away either; in sum, any additional wisdom regarding the <0DPRM would be much appreciated. 
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